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Light is produced when a scanning tunneling microscope is used to probe a metal surface. Recent 
experiments on cobalt utilizing a tungsten tip found that the light is circularly polarized; the sense 
of circular polarization depends on the direction of the sample magnetization, and the degree of 
polarization is of order 10 %. This raises the possibility of constructing a magnetic microscope with 
very good spatial resolution. We present a theory of this effect for iron and cobalt and find a degree 
of polarization of order 0.1 %. This is in disagreement with the experiments on cobalt as well as 
previous theoretical work which found order of magnitude agreement with the experimental results. 
However, a recent experiment on iron showed 0.0 ±2 %. We predict that the use of a silver tip 
would increase the degree of circular polarization for a range of photon energies. 

PACS numbers; 61.16.Ch, 78.20.Ls, 73.20.Mf 



I. INTRODUCTION 



The last ten years have seen a rapid development of 
the field of magnetic dichroism, especially where the re- 
sponse of a system to left and right circulajdy polarized 
light is probed; so called circular dichroism. El For a mag- 
netic material this is magnetic circular ^chroism (MCD). 
Experiments involve x-ray absopptionjj as well as stan- 
dard photo-emission techniques. 

The possibility of a new method of microscopic mea- 
surements of surface magnetism has been suggested by a 
recent experiment .□ In this experiment, circularly polar- 
ized light emitted from a scanning tunneling microscope 
(STM) was observed when the surface of a ferromagnetic 
material (Co) was probed with a W tip in a longitudi- 
nal configuration (see Fig. |l|, the applied magnetic field 
is parallel to both the surface plane and to the plane of 
light detection) . The handedness of the circular polariza- 
tion was found to depend on the direction of the applied 
field and the degree of polarization was between 5 and 
10%. The results of Vazques de Parga and Alvaradocl 
looks, at first sight, to have been corroborated by the 
theoretical work of Majlis et al.l3 In principle such an ef- 
fect, if confirmed, should make it possible to map the 
magnetic microstructure of a surface by measuring the 
circular polarization of the emitted light while scanning 
the surface because a STM tip provides very good spa- 
tial resolutionJj However, _a more recent experiment of 
this type, by Pierce et al.El carried out on Fe with a W 
tip, found 0.0 ±2 % circular polarization. 

In this paper we calculate the degree of circular polar- 
ization for Fe and Co. We find two contributions to the 



circular polarization. The first is due to the Kerr rotation 
of the light emitted in the tunneling process. The second 
contribution is due to the polarization of the scanning tip 
by the electric field of the emitted light. The polarized 
tip radiates and the radiation undergoes a Kerr rotation. 
This second contribution depends on the polarizability of 
the tip as well as the dielectric properties of the sample 
and can be significantly larger than the first contribution. 

We obtain results which are between one and two 
orders of magnitude smaller tiian that measured by 
Vazques de Parga and Alvaradjoa but consistent with the 
measurements of Pierce et al.lD. Both workers used W 
tips, but we find that for a Ag tip, Co and Fe produce a 
larger degree of circular polarization for a range of pho- 
ton energies. 

There are other sources of circular polarization that 
are not magnetic in origin. For example, experiments by 
Vazques de Parga and AlvaradoB and theory by Anisi- 
movas and Johanssono show that an asymmetric STM 
tip can produce a degree of circular polarization on the 
order of 10 %. 

In section || we develop a theory for the degree of circu- 
lar polarization produced by an STM tip in the presence 
of a magnetic sample based on the magneto-optic Kerr 
effect. In order to understand the physics of the mag- 
netic circular dichroism in these experiments the tip is 
modeled by a dipole. The purpose of the dipole model is 
solely to help to understand the physics and is in troduced 
mainly for pedagogical reasons. In section III the theo- 
retical description of the tip is improved for the purpose 
of obtaining reliable numerical estimates of the degree 
of circular polarization. These results are discussed in 
section IV. 
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II. THEORY 

In this section we obtain an expression for the circular 
polarization of light emitted when an STM tip scans a 
magnetic material. The calculation is divided into five 
parts: (A) We first describe the experiment and express 
the results in terms of Stokes parameters. (B) The elec- 
tromagnetic fields in the tip region are related to the 
fields at the detector. (C) A model for the tip is intro- 
duced which makes the problem tractable. (D) The local 
field outside the tip is related to the tip polarization, in 
the presence of a magnetic substrate. (E) The different 
contributions are assembled and we obtain a theoretical 
expression for the changes in the circular polarization due 
to a change in the direction of the applied magnetic field. 



A. Theoretical expression for MCD 

In this section we relate the MCD to the field ampli- 
tudes at the detectos. The experiments by Vazques de 
Parga and Alvaradoo and Pierce et al.0 used the longi- 
tudinal configuration, i.e. the applied magnetic field is 
parallel to the plane of incidence and in the surface plane 
of Co(OOOl) [thin film grown on Au(lll)] and Fe(OOl) 
(whisker), respectively. Light was detected at an an- 
gle of 30° measured from the surface and tungsten tips 
were used in both experiments. The emitted radiation 
showed circular polarization which changed when the ap- 
plied magnetic fieldrnwas reversed. Whereas Vazques de 
Parga and Alvaradoo used a fixed quarter wave plate and 
cawry out the analysis by a linear polarizer. Pierce et 
al.Q use a rotating quarter wave plate and a fixed lin- 
ear analyzer. In both experiments control measurements 
were performed on clean Au(lll) samples and yielded 
no change of the polarization of the emitted light upon 
reversal of the external magnetic field. The background 
due to geometric details of tip-sample junctions such as 
those discussed in Refs. ^ and ^ and any residual dichro- 
ism of the view port of the UHV system were removed 
by reversing the magnetization of the sample. 

The results of the two groups were very different. 
Pierce et aLQ found no magnetization dependent circu- 
lar polarization within an experimental uncertainty of 
±2 % whereas Vazques de Parga and AlvaradoQ found 
values of the order 5-10 %. An earlier experiment us- 
ing a Ni tip, thus injecting spin-polarized electrons, and 
a Ni polycrystal sample also showed a large MCD upon 
reversal of the magnetization of the tip.E3 

Expressing the degree of polarization of the emitted 
light in terms of Stoke's parameters we have 



P -P 



So 



(1) 



The superscripts on S3 and p indicate whether the mag- 
netization is parallel or antiparallel (— ) to the sur- 
face projection k, of the photon wave vector. Following 



the definition given by JacksonS and concentrating on 
p"*", left-circular polarized light (positive helicity) has the 
polarization vector 6+ = (p — is)/\/2, whereas for right- 
circular polarization e_ = {p + is)/\/2, where, see Fig. |^, 
s = X and p = 9 are s- and p-polarization unit vectors. 
In terms of the electric field ET„ at the detector caused 
by the source at the tip (T), the Stokes parameters are 
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where the last approximation follows because |i?J| is nor- 
mally much smaller than for the particular set-up 
we are considering. Hence, the experimentally measured 
quantity is 



2Im 



i?J (detector) 
Ej (detector) 



(M^-M). (3) 



The basic physics behind these equations is that tunnel- 
ing electrons undergoing inelastic events via spontaneous 
emission produce a radiating field in the vicinity of the 
tip. This leads primarily to emission of p-polarized light, 
but when reflected in the surface of the magnetic sample 
it also gives rise to a small field component in a direction 
parallel to the surface due to the Kerr effect. We have 
thus expressed the measured MCD in terms of the field 
amplitudes at the detector. To facilitate the ensuing cal- 
culations, we will next relate these tip-induced fields in 
the detector region to complementary fields generated by 
sources at the detector. 



B. Reciprocity theorem 

In this subsection we use the reciprocity theorem to 
reformulate Eq. (|^). The electrons tunneling inelasti- 
cally between tip and sample can emit photons and are a 
source of electromagnetic radiation. This radiation cou- 
ples to tip and sample and is fioaJly detected far from the 
tip. We have previously founcO that it is convenient ta 
use the reciprocity theorem of classical electrodynamicsll^l 
in such a situation because it allows the radiated field 
to be approximately determined by a nonretarded cal- 
culation if the wavelength of the emitted light is large 
compared to the relevant tip extension. This theorem 
essentially states that the result of a measurement is un- 
changed if the source and field points are interchanged. 
Here, the reciprocity theorem can be written as 
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J dVEj(^) Jf (x) = ^ I dVEf{^) Jf{^), (4) 



where j denotes the components x, y, and z. The current 
at the tip (detector) is the source for the electro- 
magnetic field E'^(-°) . 

Equation (Q) is valid for media with time-reversal sym- 
metry. However, a magnetic material does not fulfill 
this condition and one has to use a modified reciprocity 
theoremtj whereby Eq. (^) can still be used provided the 
true medium is replaced by its complementary medium 
(one with reversed magnetic field). If the dielectric ten- 
sor of the real medium is eij, then the complementary 
medium has a dielectric tensor ej- which is the transpose 



of Cy-; 



The dielectric matrix has the form 



65 (w) eicos0 — eisin7sin( 

— eicos0 es(t^) eiC0S7sin(/i 

eisin7sin(/) — ei cos7sin(/) es(w) 



(5) 



where the notation ei((j-') = iQes{^) is sometimes used 
with Q being the so called magneto-optical constant and 
£5(0;) is the substrate dielectric function. The angles 
and 7 specify the direction of the applied magnetic field 
with respect to the surface normal and the plane of in- 
cidence. We see that the complementary medium corre- 
sponds to changing the sign of the off-diagonal compo- 
nents. 

To a high degree of accuracy, the current between the 
tip and sample is spatially well-localized and perpendic- 
ular to the substrate. Thus we write 



Jj(x) = zJo{p,z) 



(6) 



where z is normal to the surface, pointing inwards. The 
relative independence of on position in the surface-tip 
region has been verified by numerical calculationcj and 
is also a key feature of the models used in this paper. 
We can take out an average value of in the region 
between tip and sample and write the right hand side of 
Eq. (I) as 



E^ j dVJo{,p,z) 
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E^ is the perpendicular component of E^, with respect 
to the surface plane, i.e. in the z-direction. The elec- 
tronic current from the STM flows primarily in the nor- 
mal direction with respect to the surface, thus only a 
perpendicular component of a "detector-generated" field 
can couple to that tunneling current. 

On the left hand side of Eq. (^, we insert a current 
source with two components, 



J^(x) = nj„(5(x - X£,), 



(8) 



where x/j is the detector position and fi corresponds to 
either the direction s or p. Notice that these currents 



generate two different ijf', which we write as -E" in what 
follows (n=s or p). Eqs. (|), (0) and (|) yield 



(detector) j„ ^ joE^ 
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Combining Eq. (|^) with Eq. (|^) and also including the 
complementary-medium sign change we obtain 
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and p~ is p^{M —M). In the last line we have re- 
placed js^p by the corresponding incoming field strengths 
Ef^p since s-polarization and p-polarization represent two 
orthogonal polarization states. Thus the reciprocity the- 
orem makes it possible to express the tip-generated field 
at the detector position (i?J (detector)) as the fields in 
the tip region generated by incoming s- and p-polarized 
waves. 



C. Determination of fields at STM tip 

In this section we model the STM tip in such a way 
as to include the main physical effects and to allow the 
development of a formalism for determining the fields at 
the tip. The system of tip and sample is a difficult one 
to treat for several reasons. Even if the tip was perfect, 
in the sense of having a well-characterized geometrical 
shape, the resulting electromagnetic field problem has 
relatively low symmetry. Consequently, we model the tip 
as a polarizable sphere of radius R, with a scalar polar- 
izability ao(^)j situated a distance d ~ R + D from a 
surface. Then we replace the sphere with a dipole at d. 
In this way we include the relevant physics, such as a rela- 
tively constant field in the region between tip and sample, 
while making the problem tractable. Since the circular 
polarization is a ratio between two quantities, we hope 
that a simpler model can capture the main features of 
the emitted light. A dipole model waa^also used to study 
the polar Kerr effect by Kosobukinlj in the context of 
near field optics. In section III we improve upon this by 



using the full sphere in the calculations. We consider the 
non-magnetic situation in this s ection and introduce the 



magnetic substrate in Sec. [I D 



Consider a sphere centered at d=(0,0,-d) where -d is 
the position of the sphere outside a metal surface whose 
optical refiection can be described in terms of its Fres- 
nel refiection coefficients ps and pp for s-polarized and 
p-polarized light, respectively. The tip (sphere) above 
the surface is replaced by a point polarizable dipole with 
polarizability ao{uj). The total electromagnetic field at 
the dipole position E(d,w) can be divided into two parts 
E'^''* and E*: E^'^' is the solution to Maxwell equations 
with an incident electromagnetic field and no tip present 
while E* is a solution when we have no incoming elec- 
tromagnetic field but the induced field at the tip plays 
the role of a source term. Assume the point dipole has 
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an induced dipole moment P=(P||,Pl). The solution to 
the E' problem can be simplified if we decompose the 
induced field at the tip in Fourier components parallel to 
the surface (k) and note that they play the role of in- 
coming electromagnetic fields analogous to the situation 
in the E^'^^-problem, (however, we have to sum over all 
possible parallel wave vector components to get the to- 
tal field). The (near) field from a dipole can be Fourier 
decomposed according to 



E(x,c^) = 



(2^ 



E(k, z, Lu)e^ 



(11) 



where E(k,z,w) is the analogue to an "incoming" elec- 
tromagnetic field. In our case the lowest order magnetic 
component is smaller than the electric field by a factpc 
Lud/c ^ 1. The total field from the dipole and its imageO 
is expressed in components parallel and perpendicular to 
the surface as 



and 



Si(k,z,cj) = -zk-E||(e'P^ -Hppe-'P^) 



(12) 



(13) 



where ps and pp are the reflection coefficients for s- and 
p-polarized light scattered from the surface.. In Eqs. (|l^ 
and (|l|) we have introduced 



E* 



— e'P''[-kP^+pP|| + — TP|| 
p p 



(14) 



and the transverse projection operator in the surface 
plane is T = 1 — kk, where k is a unit vector along 
k. Furthermore -I- fc^ = 9^ = /<? for the wave vector 
q — (k,p) of the incoming field. In order to use the reci- 
procity theorem above we expose the surface and the tip 
to an incoming electromagnetic field E'"'^ . In the absence 
of the tip the total field would be E'''^* = E'"'^ + 
where E''"'' is the reflected field. Upon introducing the 
dipole, it will develop an induced polarization 



P(d) =ao(c^)[E°'^*(d)+E*(d)] 



(15) 



where d =(0,0, -d) is the position of the dipole and E*(d) 
is the image field of the pularized tip due to the presence 
of P. One can show thattj (see also Appendix ^): 



E*(d) = i^||P|| +^_lP_lz 



(16) 



for P in the surface plane (P||) or perpendicular to it 
(Pj_). We have defined the feed-back, or image functions. 



dkk- 



-ip 



[q^Ps{k, uj) - p^Pp{k, uj)] (17) 



and 



dkk^ 



-ip 



-[pp(fc,w)] 



(18) 



after performing an angular integration in the surface 
plane. The non-retarded limit for F|| and Fj^ is obtained 
by letting c — > cxd with the result that p is replaced by 
ik. 

Equation (15) is a self-consistency condition on the in- 
duced dipole moment representing the tip. Solving for P 
we obtain 



1 - ao{uo)F, 



.i?r'(d), 



(19) 



where i=|| and _L. The denominator in Eq. ( [L9[ ) can be 
included with the field to form an effective field acting 
on the unperturbed tip (dipole) or it can be included 
with the bare polarizability ao to form an effective po- 
larizability. For special frequencies, the coupled system 
exhibits resonances when Re(aoFi)=l and Pi canJae very 
large. An explicit demonstration of this is foundEII for a 
sphere outside a surface. Combining Eqs. (15) and ( p^ ) 
we obtain the total field at the dipole position. 



E*°'(d) - P(d)/a„ 

= E'='^*(d) +E*(d) 

= E'=''*(d) -FG||Ef*(d) 



(20) 



where G\\^^_W^^ is the field at the tip, due to the image of 
the tip produced by P\\,±- We have introduced an image 
factor Gi defined as 



ao(cj)F||, 



1 - ao{uj)F\i 



(21) 



For a substrate characterized by a frequency dependent 
dielectric function, 65(0;), we find in the non-retarded 
limit: 



Fi = 2F« = 



1 es(^)-l 
4d3 es(t^) -I- 1' 



(22) 



In Eq. (I21I), Uo contains information about the dipole 
resonances and F|| _l contains information about the sur- 
face (sample) resonances (£s-|-l=0 corresponds to surface 
plasmons). From Eq. ( pl| ) we see how the tip and sam- 
ple couple to yield new eigenmodes at the poles of G||_^, 
and also a possible field enhancement. In this respect our 
model contains all the features ot-more refined calcula- 
tions for the STM configuration.EJ For the dipole (tip), 
we use the polarizability for a sphere of radius R and 
dielectric function eT{uj): 



ao{io) = R 



1 



eT(w) +2 



(23) 



In this way one can also include a more accurate dielec- 
tric response of the tip, e.g. using measured values for 
et- The factor 2 in the denominator is actually (-1+1/n) 
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where n is the so called depolarization factor. Thus we 
could also mimic different tip shapes by choosing differ- 
ent values of n. 

In the following discussion we retain Eq. (^ ) as a 
generic form for the resulting field at the tip position 
when an incoming field is incident on the tip and the 
sample. Notice also that the only boundary condition 
matching is done in the absence of the tip. After this the 
self-consistency condition for the induced dipole moment 
[Eq. (|l5|)] adjusts the total field strength appropriately. 



D. Local field near tip with magnetic substrate 



In section II B we related the fields at the detector to 



the fields of the inverse problem where light is scattered 
from the tip. In the previous section we showed how the 
incoming field is affected by the presence of a tip out- 
side a non-magnetic solid. Furthermore, we expressed 
the fields in terms of the total field outside the substrate 
in the absence of the tip: E™*(d) in Eq. (pO|). Now we 
will calculate E°''*(d), and also consider the changes in 
the reflected dipole field, in the presence of a magnetic 
sample. We first address the field felt by a tunneling 
electron that undergoes an inelastic event leading to the 
light emission. 

In the previous section we calculated the induced po- 
larization at the tip in a self-consistent manner and ob- 
tained the field at the tip position itself. A tunneling 
electron will not only feel the incoming field and the re- 
flected field of the incoming light, but the tip-induced 
image field as well as the direct dipole field from the tip. 
In the region between the actual tip and substrate these 
fields vary only slightly since the dipole is far away from 
the surface, R>>D. We can write the field acting on a 
tunneling electron as (i=||,_L): 



E, = FoPi = Foao(l + Gi)(l + pi)E!" 



(24) 



where all fields represent average values in between tip 
and sample. Fo is a the dipole factor which gives the 
field strength for a given polarizability P. Notice that 
if we have no surface present E — Fo,dircctaoE'"'^, where 
Fo, direct IS the direct part of Fq. We see that (1+G) 
plays the role of a field enhancing factor for the incom- 
ing and scattered electromagnetic fields from the surface 
(I + p)E™'^. cto is a property of the tip, G depends on 
both the sample and the tip, and p is a property of the 
sample. 

To determine P we first have to calculate the reflected 
field from a magnetic surface and such a reflection in- 
volves the magneto-optic Kerr effect. For general angles 
between applied magnetic field, surface plane and plane 
of incidence there is an excellent treatment by Zak et al.E£l 
which is very useful when dealing with light reflection 
from a magnetic solid. For a general magnetic configu- 
ration the proper form of the dielectric matrix was given 
in Eq. (H) . The off-diagonal elements of this tensor carry 



the information about the Kerr effect. The origin of the 
non-diagonal components is a coupling between the spin 
of the electrons in the solid and their orbital momentum 
due to the atomic potentials (spin-orbit coupling). The 
first theoretical calculationa.of this effect were carried out 
by Hulmellj and Argyres.Ej Equation (H) assumes that 
the dielectric tensor is diagonal when ei(ti^) — (this ap- 
proximation can be relaxed but it does not influence our 
final conclusions). 

For the applied magnetic field in the surface plane we 
can use Eq. (|^) for the general dielectric tensor with 
= 7r/2 and one angle (7) suffices to specify the di- 
rection of M with respect to the plane of incidence. The 
dielectric matrix in Eq. (^ then simplifies to 

es(w) — ei(ci;) sin7 \ 

es(^) ei(a;)cos7 . (25) 

£1(0;) sin 7 —61(0;) cos 7 es(ti^) / 

We use this result in the Zak matrix multiplication 
method, together with an expansion to first order in 
£1(0;), since the off-diagonal elements are small com- 
pared to es(tf) (I t\{iS) I ts{^) 1^1)- However, these ma- 
trix elements provide a coupling between s-polarization 
and p-polarization leading to a non-zero Kerr rotation. 
For s-polarized incident light (along x-direction) , i?™'^, 
our analysis gives the reflection coefficient for s-polarized 
light to lowest order in as 



V~Vs 



(26) 



P + Ps 

which is the standard Fresnel result .0 — fc^ and 

for q > k, k — qsind with 6 being the angle of inci- 
dence. For k >> q, p ik. ps — \/ q'^es{(jj) — fc^, with 
es(w) deflned above. Finally q=w/c, where c is the light 
velocity. Apart from the direct reflection of s-polarized 
light there is a small but crucial conversion from s- to 
p-polarized light with a reflection coefficient 



Pps = 



kpq 



Ps{p + Ps){esP + Ps) 



sin 7 -I- 0{el) = pL sin 7 
(27) 



It is proportional to £1(0;) since it comes from the off- 
diagonal response. For large es, Pps oc ei/el- With Zak's 
field conventions, Eq. ( p7| ) corresponds to a reflected fleld 
amplitude, — p^s^'"'^ sin^, in the z-direction. It is clear 
from Eq. ( pTj) that a general direction of the applied mag- 
netic field in the surface plane corresponds to replacing 
£1(0;) in the strictly longitudinal configuration (7 = 7r/2) 
by ei(a;) sin 7 for the s-polarized case. Changing the di- 
rection of M (7 — > 7 -I- tt) changes the sign of pps as it 
should. 

Repeating the same analysis as above for p-polarized 
incoming light we find, to lowest order in the off-diagonal 
elements: 



Pp 



<^sP - PS 
<^sP + PS 



+ Pt cos 7 4- 0{e\) 



(28) 
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with 



Pt 



—2ti{uj)pk 
{es{uj)p + Ps) 



(29) 



The corresponding reflection coefficient for p- to s- 
polarized conversion is psp = —Pps- The first term in 
Eq. (|2^ ) is the standajaJ Fresnel reflection coefficient for 
a non-magnetic solid.Ell For later use we need the non- 
retarded limit of equations (28) and (E9), viz.: 



Pn = 



esjuj) - 1 
es(w) + l 



where 



Pt 



p°rp cos7 + 0(ei) 



2iei(aj) 



(30) 



(31) 



The first term in Eq. (|30| ) is the classical image factor for 
a solid with dielectric function es(w). p-. 

Making use of the formalism of Zak et al.EJ, we find 
that for an irradiated magnetic surface we can make the 
following replacements with respect to the non-magnetic 
situation: 



Ps 



PL sm 7 



and 



Pp^ Pp + Pps ^ Pp + PL sin7 + Pt cos-y. 



(32) 



(33) 



PL is defined in Eq. (|27|) and px is defined in Eq. (^9|). 

The change in reflection factors for a magnetic surface 
compared to the non-magnetic situation will also affect 
the reflected field from the tip. We are interested in the 
near-field and therefore take the non- retarded limit of p^ ; 

pV- 



PL 



qei{uj) 



2fc(es(w) + 1) 



(34) 



to lowest order. Notice that pp has a finite value [p° + 
/9^cos7, cf. Eq. (|30|)] in this hmit and that ps vanishes 
as 0((q/k)^). In what follows, we neglect pi oc q/k com- 
pared to p^ since fc 3> w/c in the non- retarded limit. Re- 
peating the previous treatment for a magnetic material 
one finds [by letting Ps and pp be transformed according 
to Eqs. (|2|) and (^)] that Eq. (^ is replaced by (see 
Appendix^): 



E|r(d) = F||P|| - FTP±{i X M) 



(35) 



and 



E'°''{d) = + FtP\i • (z X M) (36) 



where M is a unit vector in the surface plane in the di- 
rection of M and z is normal to the metal and directed 
towards it. The coupling factor due to the off-diagonal 



response of the medium, Ft, is given by (non-retarded 
limit) : 



Ft=y^ J dkk'^p?p(k,Lo)e- 



■2kd 



1 ei(^) 
4d3 {es{Lo) + If 



(37) 



using from Eq. (|3l| ) in the last line. The physics be- 
hind the structure of the above equations is the following. 
A perpendicular dipole [Eq. (|35|)] provides an electro- 
magnetic field which is reflected in the surface and gives 
an induced electric field and dipole component parallel 
to the surface (due to Kerr response, through Ft) and 
perpendicular to the magnetization M. A parallel dipole 
likewise is reflected and provides an induced perpendic- 
ular field and dipole, due to Ft- If both M and P|| are 
parallel there is however no such contribution. 



E. Theoretical results 

We have now developed all the necessary ingredients 
for calculating the circular polarization from Eq. (p^). 
First consider the case of incident s-polarized light at an 
angle 9. A field Ef'^s (along the x-direction, s = x) sets 
up a parallel polarization of the tip 

P|| =ao[(l + Ps)Sri + F||P|| -FtP±(zxM)] (38) 



Eq. (y_5|) and Eq. ( pS]) have been used to derive Eq. (|3q). 
There is also a perpendicular polarization 

Pi_ = ao[- sin0p(,;*^) ET^ + F^Pi_ + FtP\\ ■ (z x M)]. 

(39) 

In the first term we have indicated that the conversion 
coefficient pps is now, due to the exact form of the reci- 
procity theorem, to be evaluated in a situation which is 
the same as the one considered above if we change the 
sign of M. The Kerr coupling between the induced dipole 
moments of the tip is included through the coupling func- 
tion Ft- In Eq. ( |39| ) the first term is the z-component 
Kerr field set up by the incoming s-polarized light. The 
second term is the image from the tip-induced polariza- 
tion perpendicular to the surface. Finally the third term 
is the image from the parallel induced dipole set up by the 
incoming s-polarized wave [Eq. (^8|)] which is converted 
to a perpendicular component by the non-diagonal re- 
sponse of the substrate. The latter is described through 
Ft which is given above. Neglecting Ft in equation (|3^), 
substituting the resulting expression in equation Jj33) for 
Pji into equation (|3^) and using equations (|l5|), (|20| ) and 
( pl| ) gives the field experienced by an electron between 
tip and sample to lowest order in Ft [cf. Eq. (|2^)]: 

Ft = {l + G^)iaoFo) 

x[pL sine - {I + p,)G\\K{lu)]EI'"' sinj (40) 



due to the incoming field E].' 
— sin 7 and defined: 



We have used x- (z x M) 



6 



K{uj) = FtIF\\ = 



2ei(w) 



1) 



(41) 



Note that the right hand side of equation ( pO| ) is propor- 
tional to sin 7, the orientation of the magnetization in 
the surface plane; longitudinal polarization providing for 
the maximum field strength. 

El in equation ( ^0| ) is the major part of the average 
field in the z direction in the narrow region between tip 
and sample, created by the incoming s-polarized field of 
magnitude i?™^. The first factor in Eq. ( ^0|) is under- 
stood as follows; an incoming s-polarized field of magni- 
tude EY"^ undergoes a Kerr rotation and a field propor- 
tional to PL is created in the z direction. The tip acquires 
a z-component of polarization and radiates. This radia- 
tion is also reflected back to the tip by the surface, thus 
the tip sees its own image giving a contribution, G±pL- 
The total field at the tip is the original field plus the im- 
age field. Thus, G± is an image (enhancement) factor 
due to the polarization of the tip in the z direction. The 
quantity G± depends on the dielectric properties of the 
tip, the metal, the distance between them, and on the 
geometry of the tip. 

The second term in equation ( ^ ) is explained as fol- 
lows; (1 + ps)E™'^ is the field at the tip due to light that 
falls directly on the tip plus light that is refiected from 
the metal surface. The tip is then polarized in a direction 
parallel to the surface and it radiates. This s-polarized 
radiation is reflected in the surface and undergoes a Kerr 
rotation so that it develops a z-component. The term 
G\\K is the analogue of pL in the first term. Once a field 
is created in the z direction it is enhanced by the factor 
(1 + G_l) in front. 

Performing the above calculation for an incoming p- 
polarized wave with amplitude i?™"^ and working only to 
zeroth order in the off-diagonal dielectric matrix (since 
both terms in Eq. ( ^0[ ) above for s-polarized light are 
already of first order in ei(w)) we find: 



-(1 + G^){aoFo){l + pp) sinOE';' 



(42) 



where snvOE™'^ is the field that falls directly on the tip 
and Pp sin OE^'^ is the field at the tip that is reflected 
from the surface. The total fleld is enhanced by the same 
factor (1 -I- G±) as discussed above. 

With the use of Eqs. (|), (|l^, (|40|), (||) and (||) this 
leads to the following expression for the magnetic circular 
dichroism, to lowest order in ei{i-u): 



p — p =4 sin 7 Im 



PL 



Gil 



Pp 



Pp 



K{lo) 



(43) 



Note that the light enhancement factor G± from the per- 
pendicular fleld component drops completely out of the 
problem as does the dipolar factor Eq so the ratio does 
not depend on where i n th e junction the light emission 
takes place. In section III we will make a more realistic 



estimate of Eq. ( p3| ) for the magnitude of the magnetic 
circular dichroism. However, let us already here say that 
the fact that G±^ drops out of the dipole-model calcu- 
lation gives one hint to why this model, as we will see, 
yields results for the magnetic circular dichroism that 
are of the right order of magnitude although the field en- 
hancement (described by G±^) in the dipolar model may 
be very different from that calculated in the improved 
model. 



III. IMPROVED THEORETICAL MODEL 

Here we outline a calculation of the tip-induced MCD 
signal within a model geometry where the STM tip is 
represented by a sphere characterized by a bulk dielectric 
function 67^(0;). This allows for a much better description 
of the tip polarization. Most of the calculational details 
are deferred to Appendix 

We set out to determine E^ in Eq. (|l0|), and this is 
achieved via the following four steps, (i) The incoming 
s-polarized wave is reflected by the sample surface. This 
yields a total fleld E'^'^'^ parallel to the surface, (ii) That 
fleld in turn drives the model tip so that it sends out a 
fleld that is reflected back and forth between the tip and 
sample. In this step, we need to extract the part of the 
fleld that the tip sends onto the sample, (iii) Next, due to 
the off-diagonal elements of the sample dielectric tensor, 
part of the electric fleld parallel to, and incident on the 
sample is converted to a field perpendicular to the surface 
as follows from Eq. (pO|). It is only at this stage that the 
magnetic properties of the sample enters the calculation, 
(iv) In the last step, we calculate the degree to which the 
converted electric fleld is enhanced inside the tip-sample 
cavity. 

Step (i): With an s-polarized wave incident from the 
right (positive y) with electric field E™"^, the Fresnel for- 
mulae yields a total field just outside the sample surface 
given by 



Er(l + p,)=fi?i' 



2p 



p + ps 



(44) 



Thus before introducing the model tip into the problem, 
we have an electric field outside the sample that can be 
described (in the non-retarded limit) by the scalar po- 
tential 



= -xE° 



(45) 



where E'=''^ = E'f'''[2p/ (p + ps)]. 

Step (ii): Once the model tip is introduced into the 
problem clf^^ alone is no longer a solution of Laplace's 
equation in the region above the sample, instead another 
contribution </)ind has to be added. Using the appropri- 
ate boundary conditions for the E and D fields at the 
sample and tip surfaces 0ind can be determined. In the 
following, we only want to keep the part of (^ind that the 
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f sends onto the sample. As we will see in Appendix 
this separation can be done by a simple inspection of 
the solution. 

Step (in): We proceed to find the field that is reflected 
from the sample surface due to the second term in En . 
®; this is the converted field E™'^^. Equation (|g) 
defines a surface response function x(k,a;) — p^, which 
in terms of incident and reflected electrostatic potentials 
is defined as the ratio [—(fy"^^i\<i,Lo)/(f)™'^{}<i^uj)\. A fur- 
ther analysis shows that within a non-retarded treatment 



there is a local relation between x ■ EJJJ|j and z ■ E™"^, 



2ei(c^) 



(46) 



Step ( iv): In this final step, we calculate the enhancement 
of the converted field due to the presence of the model 
tip. The converted field discussed above can be repre- 
sented in terms of a scalar potential (ff . Again, with the 
tip present, (jf alone does not solve Laplace's equation; 
it must be supplanted by another contribution (jf^^^- The 
calculation determining (/if^^ is completely analogous to 
the one carried out in Ref. |l2|, the only difference being 
that (jf is the driving "force" in the present case. Having 
found (jf and (t>'{^^, we evaluate the corresponding elec- 
tric field on the symmetry axis. This is the tip-induced 
contribution to El appearing in Eq. (p^). 



IV. NUMERICAL RESULTS AND DISCUSSION 

We now present numerical results for the dipole model 
of the tip based on the expression in Eq. (|43|) , and for the 
sphere model of the tip discussed in Sec. |II4 We use ex- 
perimental optical data for the dielectric functions of the 
tip and the sample. The off diagonal matrix elements of 
the sample dielectric function, ei(w), are obtained from 
magneto-optic Kerr effect measurements. The literature 
contains many detailed calculations and measurements 
of the Kerr effect (no tip present). The Kerr effect is an 
optics effect caused by the spin-orbit interaction that was 
discovered in the last century. The spin-orbit interaction 
is small in Fe and Co because the orbital momentum in 
3d-metals is small. It is only in. this century that a mi- 
croscopic theory has emerg ed.E3T2l The Kerr effect has 
recently been calculated by a number of groups for a va- 
riety of elements,aad compounds (see, e.g. Gasche et al. 
and Delin et al.E3'Eil). Similarly, on the experimental side 
there have been a number of measurements from those of 
Krinchik and Artem'evE3, whose results wo-use to obtain 
ei((jj), to the recent results of Weller et al.E3 

Equation (|3|), derived for the dipole model, consists 
of two terms; the first corresponds to the direct Kerr ro- 
tation by the sample and the second to the Kerr rotation 
of the light produced by the radiating polarization of the 
STM tip as previously discussed. We will refer to the 
two effects as the substrate Kerr effect and the tip Kerr 
effect. 



In order to use the dipole model to make an estimate 
of the tip-induced Kerr effect we calculate G|| from Eqs. 
(^), (p2[), and (23). The distance between the surface 



and the tip (i.e. from the surface to the sphere) is very 
small cornpared to the tip radius so that d ^ R in Eqs. 
( p^ ) and ( |2^ ) and G|| is determined by using experimental 
values for the dielectric functions of the tip and sample. 
When comparing our results in Eq. (^ ) with the experi- 
ments of Vazques de Parga and Alvaradoo and Pierce et 
al.Ll, we immediately recognize that in the longitudinal 
configuration that they use, s and M are perpendicular 
to each other so that sin 7= I. 

Considering the improved (sphere) model, the quanti- 
ties entering Eq. ( were calculated as outlined in Sec. 
Ill and Appendix ^ As for the model geometry we here 
used a tip-sample distance D =5 A, and the tip radius R 
was set to 300 A. 

The results of the calculations for the degree of circular 
polarization, p+ — p", are presented in Fig. Figure ^ 
(a) displays the results obtained with a Co sample, while 
Fig. H (b) shows the results relevant for a Fe sample. In 
each panel, the results for the substrate Kerr contribu- 
tion (one single curve) and the tip-induced contribution 
(four curves) to the degree of polarization are presented 
separately. The tip-induced degree of polarization has 
been calculated for Ag and W tips, respectively, using 
either the dipole model or the sphere model. The sphere 
model gives a larger degree of circular polarization than 
the dipole model, because it is a better description of the 
local electromagnetic interaction between the tip and the 
sample. In our calculational schemes, the tip-induced 
contribution ultimately results from a Kerr rotation of 
the electromagnetic fields incident from the tip onto the 
sample. The strength of the incident field is determined 
by how strongly the tip is excited by the incoming s wave 
as well as by waves reflected back and forth between the 
tip and sample. The sphere model allows for a more com- 
plete treatment of the repeated reflections than the dipole 
model, and therefore gives larger values for p+ — . The 
dipole approximation works fairly well because the circu- 
lar polarization is given by the ratio of p and s electric 
fields. 

The effects of the tip-sample interaction is particu- 
larly pronounced for a silver tip. An isolated sphere 
has a dipole-plasmon resonance when the dielectric func- 
tion exioj) = —2. For silver, the dielectric function ap- 
proaches this value near 3.5 eV. Then the model tip be- 
comes highly polarizable and the feedback mechanism, 
the waves reflected between the tip and the sample, de- 
scribed above becomes even more effective. In connection 
with this the electromagnetic response functions that en- 
ter our calculations undergo large phase shifts so that 
p+ — p~ changes sign one or several times.cS The rela- 
tively small magnitude of the results for the MCD is due 
to the factor [eg + 1)^ in the denominator in Eqs. ( |3l| ) 
and (p6|). Since 65 for both Co and Fe is rather large this 
suppresses the MCD signal. From bulk arguments one 
could have assumed that the circular dichroism should 
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be proportional to ci/es- However the presence of the 
surface changes the magnitude of both the field going 
into the solid to be Kerr rotated and the resulting field 
going out again; in both cases with a factor l/(es + 1). 

It is immediately clear that within the model we have 
considered, there is no explanation for the large values 
for the degree of light polarization found by Vazques de 
Parga and Alvarado.u The calculated results are an or- 
der of magnitude or more smaller than the experimental 
results found in Ref. ^ Even though we use a rather sim- 
ple model for the geometry as compared with the com- 
plicated, and to a certain degree unknown geometry of a 
real STM tip, we cannot see how this could make up for 
the very large difference between experimental and theo- 
retical results. Varying the geometry parameters within 
reasonable limits (within a non-retarded formulation the 
results depend only on R/D) can change the calculated 
degree of polarization by a factor of 2 at most. This is 
also the case when using different sets of optical data for 
the dielectric functions entering. On the other hand, our 
calculated results are consistent with the experimental 
results found by Pierce et al.l3 on Fe samples. Our re- 
sults are also small compared to those obtained in the 
calculations of Majlis et-al.,El some of the reasons for this 
are given in a footnote. E3 

Our results are not affected by the assumption of a 
bulk sample compared to the very thin Co (lOOA) film on 
AuQ, provided it can still be described with bulk dielectric 
data. Calculations by Moog et al.E2l indicate that there 
are no fundamental changes in actual numbers for films 
ranging in thickness between 100 Angstrom and 400, ex- 
cept a slight enhancement of the Kerr parameters. 
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APPENDIX A: 

In this appendix we show how to evaluate the differ- 
ent angular integrations required in averaging the local 
dipole field over a magnetic surface. 

We first have a look at the integrations for the non- 
magnetic case. In this instance we have three types 
of integrals appearing / fi^A:(T • Ep, J d^fc(k • E|) and 

Jd^A;Ej|. They can be expressed in terms of the three 

integrals / (fkk, J d^k(k-P\\) and / d2fc(k(k • P,, )). Per- 
forming the angular integration we find that the first in- 
tegral is identically zero. Since the second integral is Pii 



dotted with the first integral it also vanishes. This leaves 
the third integral (where if is the angular variable): 



dip 
2^ 



k(k.P||) 



(Al) 



This result is obtained by expressing k as x cos (f+y sin (p 
and k • P|| = Px cos (p + Py sin ip. 

In the magnetic case we have two extra angular fac- 
tors in the integrand compared to the non-magnetic 
case; siwy and cos 7. Expressed in terms of k they 
are k • M and z • (M x k), where M is a unit vec- 
tor in the direction of the applied magnetic field. Re- 
peating the same steps as above we find that all inte- 
grals can be expressed in terms of the following four 
integrals: J d'^k{z ■ (M x k)), J d'^kk{z ■ (M x k))), 
/ d2fc[(k.P||)(z-(Mxk))] and,/ d2/!c[k(k.Pj|)(z-(]VIxk))]. 
The first and last integrals vanishes since they are odd 
in k. The second integral becomes ^(z x M) using Eq. 
(Al). For the third integral we obtain: 



zvr 



1 



P||)(i-(Mxk)) = -P,|.(zxM) (A2) 



again expressing the different vectors in x and y compo- 
nents and using Eq. (Al). With the use of these results it 
is a straightforward manipulation to arrive at Eqs. ( |35| ) 
and M). 



APPENDIX B: 



This Appendix explains the calculations outlined in 
Sec. [II in more detaxL To carry them out we use 
bispherical coordinates, Ell furthermore, to facilitate the 



connection with similar, earlier calculations by usIiS and 
others,E3 we also introduce another Cartesian coordinate 
system {x',y',z') in which x' and z' are reversed com- 
pared with X and z, see Fig. ^. The bispherical coordi- 
nates {(3,a,Lp') are defined by 



a sinh (3 



x' 4- iy' 



asm a e 



(Bl) 



cosh /3 — cos a cosh /? — cos ol 

Both the sample surface (z' = 0, /? = 0) and the sphere 

/3 = /5o = ln(l + (D + a)/i?), 

are constant /3 surfaces. The length scale is set by the 
parameter 

a = \Jd'^ + 2RD. 



Let us embark on the calculations described in Sec. [II]. 
Step (i) is a straightforward application of the Fresnel for- 
mulae that yields an external potential 0'^''* = —xE'^^^ = 
j^'jTjcxt jggg (^5|)] describing the electric field of an s 
polarized wave refiected off the sample surface. In the 
bispherical coordinates this potential can be written 
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_^cxt ^ ^cosh/3 cos a 

oo 

c^e-("+i/2)l^lp„i(cosa) cos(p'. 



(B2) 



Proceeding to step (ii), we introduce the induced poten- 
tial (/)ind, and make the Ansatz 



-VS a v^cosh/3 - cos a 

OO 

;^F„(/3)P,i(cosa) cos^^'. 



(B3) 



In the sample {/3 < 0), the function i^„(/3) is given by 

K(/3) = (A„ + i?„)e("+i/2)'3, 

in the tip (/3 > /3o) 

F„(/3) = (A„e(2"+i)'5o B„)e-("+i/2)/3, (B4b) 

and finally, between the sample and tip (0 < /3 < /3o) 

F„(/3) = A„e("+i/2)^ + B„e-("+i/2)/3. (B4c) 

Thus, Fn, (j)ind, and therefore the tangential E field is 
continuous across the tip and sample interfaces. From 
the form of Eq. ( B4c ) it is clear that the field that the 
tip sends onto the sample is contained in the A„ terms; 
these fields decay exponentially as one goes away from 
the tip. 

To determine the coefficients, An and Bn we must also 
demand that the displacement field perpendicular to the 
sample and tip surfaces is continuous across these inter- 
faces. At the sample this means that 



find I 

aa 1/3=0- 



I'ind 



1/3=0+ 



(B5) 



(here £5 is the sample dielectric function), which yields 
Bn - -XsAn, (B6) 
where the sample surface response function 

es-l 



Xs 



(B7) 



Note that the contribution to D± coming from (p^'^^ is al- 
ready continuous since the external potential results from 
using the Fresnel formulae. At the tip-vacuum interface 
{(3 = Po), both (j)ind and (p°^^ must be considered in the 
boundary condition for D±^ . This yields an equation sys- 
tem involving the coefficients A„ 



(B8) 



where 



f/;= -(2n + 1) cosh/3o (e("+i/2)'3« - xs XT e-("+i/^)*) , 

+ Sinh /3o XT (e("+l/2)/3o _ g-(n+l/2)/3o j (B9a) 



= (n - 1) (e("-i/2)/3o „ ,-(«-i/2)/3„^ (B9b) 



n 

and 



{n + 2) fe("+3/2)/3o_^^^^e- 



(n+3/2)/3, 



(B9c) 



5„= -XT |e-("+i/2)'5o [sinh/3o - (2n + 1) cosh/3o] 
+ (n - i)e-("- V2)/3o _^ ^ 2)e-("+3/2)/3o | _ (^g^-^ 

One arrives at these equations through a procedure 
that is complfttftly analogous to the one used in earlier 
calculations .Il3£3 Solving Eq. (B8) for the An coefficients. 



we can determine the tip-induced potential incident on 
the sample. 



/mc 
^'ind 



-VS a E"'^^ \/cosh /3 - cos a 

00 

; ^^„e("+i/2)'3p,i(cosa)cos</7'. (BIO) 



Next, we have to consider the electric field conversion 
at the sample surface due to the off-diagonal components 
of the dielectric tensor in Eq. (p5|). In the primed coor- 
dinate system, it takes the form 

es{Lj) — ei(u;)sin7 

€s{lu) -ei(w)cos7 | . (Bll) 

ei(aj)sin7 ei(a;)cos7 es(w) 

Let us use this to calculate the modified sample surface 
response function. If a potential 



(B12) 



acts on the sample, its response yields another contribu- 
tion to the potential above the surface 



-Xs(k,w)e 



ili-p —kz' —iujt 



(B13) 



where xsO^j ^) is the surface response function. Combin- 
ing these two contribution to the potential with a solution 
to Laplace's equation inside the sample 



i\<i-p kz' —iiot 



(B14) 



and applyi ng th e usual boundary conditions, we obtain 
using Eq. (pil|) that 



Xs(t^) 



es-l 2ieisin(7 + (^'j.) 



es + 1 



(B15) 



where is the angle between k and x' . 

The second term in Eq. (B15) is our main concern; 
it governs the electric field conversion at the surface of 
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the magnetic sample. In view of Eq. (B15) we write the 
potential corresponding to the converted field as 



2iei{uj) 



sin(7 + ^'fc)Cd(k). (B16) 



Since we do not know the Fourier transform (/>!JJ|^(k) we 
cannot immediately use this relation. However, specializ- 
ing to the longitudinal configuration (7 = 7r/2) and using 
the fact that on the sample surface, — f{p) costs' {p 
is the distance to the symmetry axis), one can show that 
with M = My' there is a local relation between x' ■ EjJJJj 
and z' ■ E™"^, 



l(n-l)n 
"""""4 2n-l ' 

w„ = - 2n + 1 - ^ 

4 V 2n + 3 2n-l 

^ 1 (n+ l)(n + 2) 

^" ^ 4 2n + 3 



and 



(B21) 



'{P) 



2ei(co) 



•ES^(P) 



(B17) 



The calculation leading to Eq. (B17) also shows that 
both these fields can be written on the form /o(/o) + 
/2(p) cos 2(/3', where /o and /2 are functions of p. The 
second term obviously cannot induce an electric field 
possessing a z' component on the symmetry axis of the 
tip-sample system, so we neglect it from now on. The 
cylindrically symmetric part of the field E'^°"'^(/9) can be 
derived from a potential written as 



The final step amounts to solving for the enhancement 
of the z' component of the electric field between the tip 
and sample due to the presence of the tip. In fact, to 
obtain the final result for the degree of circular polariza- 
tion, this calculation has to be done with two different 
driving forces. On one hand, using (jf as a driving force 
we obtain to be used in Eq. ( p^ , if instead, as in Ref. 
p^ , the potential corresponding to a reflected p wave is 
used to drive the tip-sample system, we obtain E'^ to be 
used in Eq. (|l^). The two calculations can be done in 
parallel. We make an Ansatz for the potentials induced 
due to the presence of the tip in the two cases 

"thnd = %/2 a v/cosh /? - cos a ^ F^}"^^ iP)Pn (cos a) , 



(B22) 



where the superscripts c and p, respectively, indicate that 
it is either the converted field or an incoming p wave that 
plays the role of driving force. In the sample (/3 < 0), 



V2 aVcosh/3-cosa^ r„e"("+i)^P„(cos a). 


(B18) F<P\P) = A^(P)e("+l/2)(/3-/3o) _,_ ^c(p)g(„+l/2)(/3+/3o)^ 



The coefflcie nts T n in Eq. (B18) can be calculated with 
the aid of Eq. (B17). To this end we extract the cylindri- 



cally symmetric part of x' ■ EJ"|^ by taking the average of 
its values on the x' and y' axis, respectively, and rewrite 
these expressions as sums of Legendre polynomials only. 
At the expense of introducing a more complicated pref- 
actor this yields 



x' ■ E 



inc,a:'(y') 



ind 



= V8 E" 



COS a 



1 -I- cos a 



X]Q'(^'^P«(cos«), 



(B19) 



where the superscript x'{y') indicates whether the field 
is evaluated on the x' or y' axis. The coefficients 
and are given by rather lengthy expressions involv- 
ing the An coefficients. It is also possi ble t o rewrite E"^, 
resulting from the z' gradient of Eq. (B18) as a sum of 
Legendre polynomials preceded by the same prefactor as 
in Eq. (B19). From Eq. (B17) we then obtain the follow- 
ing equation system determining the r„'s 



(B23a) 

in vacuum (0 < /3 < /3o) 

Fn^'^\(3) = ylc(p)g(«-hl/2)(/3-/3o) _^^c(p)g-(«+l/2)(/3-/3o)^ 

(B23b) 

and in the tip {(3 > /?o), 

F<P\P) = (^A<P^ + B=(P)) e-("+i/2)(/3-/3o)^ (B23c) 

The external po tenti al in the case of the converted field 
is given by Eq. (BIS) 



= \/2 av/cosh /? - cos a ^ r„e"("+*)'^P„(cos a), 



whereas in the case of an incoming p wave it is 



— sm 



9 -JflE— E"'" V2a v/cosh/3-cosa 



+ WnTn+2 



where 



ei(^) 



(B20) 



esp + Pt 

00 

X ^(2n + l)e-("+i/2)/3p„(cosa), (B24) 


where the first factor oc i?™'^ is the resulting z' compo- 
nent of the E field just outside the sample surface as 
obtained from the Fresnel formulae. 
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Proceeding along the same lines as when deriving Eq. 

(B25) 



(B8), we first obtain 
and then the equation systems 
and 

C/„AP + K^Li + W'„AP+i = ^p 



(B26) 



(B27) 



determining the enhancement due to the presence of the 
tip of the converted field and the p polarized wave, re- 
spectively. The coefficients on the left hand side are the 
same in both cases, 

C/„ = (2n + 1) cosh/3o (l - XsXre-^^^^^^"") 

-XT sinh/3o (l - Xse-(2"+i)''°) ^ (B28a) 



(B28b) 



W„ - -{n + 1) (l - XsXTe-('"+3)''°) . (B28c) 
The right hand side in Eq. (|B26| ) is 

S'n = XT e-("+i/2)/3o {[sinh/3o - (2n + 1) cosh/?o]T„ 

+ne''«T„_i + {n + l)e-^«T„+i} , (B29) 

while for the case of an incident p wave we obtain 



S'P== sin^ 



2esp 



-Ef'XTe 



-(n+l/2)/3o 



(SP + Pt 

X {(2n + 1) [{2n + 1) cosh/3o - sinh/3o] 

-(2n2 - n)e'^° - {2n^ + 5n + 3)e''^°} . (B30) 

Once all the coefficients A:, and Ap have been deter- 



mined using Eqs. (B25), (B26), and (B27), we can cal- 
culate the resulting electric fields on the symmetry axis 
as 



and 



cosh f3 — cos a 



cosh /? — cos a 



dp dfi 



dp dp 



(B31) 



(B32) 



The degree of polarization of the emitted light is found 
from 



|^«-2Im 

Do 



E° 



zpout 

p 



(B33) 



so that using the reciprocity theorem, keeping in mind 
the sign change of ei discussed earlier 

p+-p-«4Im II . (B34) 
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Detector 



FIG. 1. The schematic experimental set-up in the experiments of Vazques de Parga and Alvaraddj and Pierce et al.u 
showing the relative orientation of the applied magnetic field (M) in the plane of the Co(QQ01)/Fe(001) surface (x,y plane), 
tungsten tip orientation, and optical detection axis. With the help of the so called reciprocity theorem [Eq. (^] one can relate 
the field intensity at the detector due to a current at the tip (E"^ (detector)) to the "detector-generated" field intensity between 
tip and sample (E^). The latter is easier to construct and hence makes it rather straightforward to find the field we are most 
interested in; the one in the tip region, s and p denote two orthogonal polarization directions at the detector. 
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FIG. 2. Calculated results for the circular polarization for (a) a Co sample, and (b) a Fe sample. Each of the figures present 
results obtained from both the dipole model (dotted curves) and the improved (sphere) model (full curves), moreover, results 
corresponding to Ag and W tips (as indicated next to the curves) as well as the substrate Kerr contribution arc displayed in both 
panels. In the sphere-model calculation we used R = 300 A for the tip radius and D = 5 A for the tip-sample separation. The 
sphere model in general gives a larger tip-induced contribution to the degree of polarization, nevertheless the results obtained 
with the two different methods have many qualitative features in common. 
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